arXiv:astro-ph/9706113vl 11 Jun 1997 


Stellar Populations and the White Dwarf Mass Function: 
Connections To Supernova la Luminosities 

Ted von Hippel 

Dept, of Astronomy, University of Wisconsin, Madison WI 53706 

email: ted@noao.edu 

G. D. Bothun 

Dept, of Physics, University of Oregon, Eugene OR, 97403 
email: nuts@moo2.uoregon.edu 

and 

R. A. Schommer 

Cerro Tololo InterAmerican Observatory, La Serena, Chile 
email: rschommer@noao.edu 

ABSTRACT 

We discuss the luminosity function of SNe la under the assumption that 
recent evidence for dispersion in this standard candle is related to variations 
in the white dwarf mass function (WDMF) in the host galaxies. We develop 
a simple parameterization of the WDMF as a function of age of a stellar 
population and apply this to galaxies of different morphological types. We 
show that this simplified model is consistent with the observed WDMF of 
Bergeron et al. (1992) for the solar neighborhood. Our simple models predict 
that WDMF variations can produce a range of more than l. m 8 in M#(SN la), 
which is comparable to the observed value using the data of Phillips (1993) and 
van den Bergh (1996). We also predict a galaxy type dependence of Ms(SN la) 
under standard assumptions of the star formation history in these galaxies and 
show that M#(SN la) can evolve with redshift. In principle both evolutionary 
and galaxy type corrections should be applied to recover the intrinsic range of 
M#(SN la) from the observed values. Our current inadequate knowledge of the 
star formation history of galaxies coupled with poor physical understanding of 
the SN la mechanism makes the reliable estimation of these corrections both 
difficult and controversial. The predictions of our models combined with the 
observed galaxy and redshift correlations may have the power to discriminate 
between the Chandrasekhar and the sub-Chandrasekhar progenitor scenarios for 
SNe la. 

Subject headings: galaxies: Supernovae; cosmology: distance scale 
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1. Introduction 

It is generally assumed that type la supernovae (SNe la) are thermonuclear explosions 
of degenerate white dwarfs near the Chandrasekhar limit (cf. Wheeler & Harkness 1990), 
or perhaps mergers of white dwarfs (WDs) in a binary system (Paczynski 1985). Either 
detonation or deflagration models (Arnett 1969; Nomoto et al. 1976) then produce the visible 
energy release that characterize the SN light and velocity curves. Detailed models show 
considerable differences in these scenarios (Khokhlov et al. 1993), but their large intrinsic 
luminosity coupled with the assumed universal physics involved in the Chandrasekhar mass 
limit have led to strong statements concerning the use of the magnitude at maximum 
(Ms(niax)) of SNe la as distance indicators (Branch & Tammann 1992). 

There has been some recent evidence that a moderate to large dispersion exists in 
Mg (max), however, and perhaps in the intrinsic color of these objects at maximum. This 
evidence has been revealed by the extensive efforts of several groups to obtain high quality 
and frequently sampled observations of a large number of SNe. In particular, Phillips (1993) 
has summarized high quality SN la measurements in nearby galaxies, with the result that 
the intrinsic dispersion in Mg(max) appears to be ~ 0. m 8, and correlated with the decay 
time of the light curve. Furthermore, the underluminous nature of SN1991bg (Filippenko 
et al. 1992; Leibundgut et al. 1993) is striking. This apparent dispersion in M B (max) has 
led to some discussion of different models of the origin of SN la explosions. Unfortunately, 
current measurements are insufficient to discriminate in detail between various explosion 
models (detonation versus deflagration, etc.) or origin scenarios (see Wheeler & Harkness 
1990; Iben & Tutukov 1991; Khokhlov et al. 1993). Woosley & Weaver (1994) and Livne & 
Arnett (1995), among others, have discussed explosions of sub-Chandrasekhar limit WDs. 
Pinto & Eastman (1997) examine the physics of SN la light curves in detail and use an 
analytic model to study the sensitivity of the resultant light curves to various properties 
of supernova explosions. They find that a variation in total mass can lead to a sequence 
of light curves that reproduces the luminosity - decline rate relation. Other possible 
parameters (explosion energy, 56 Ni mass, and opacity) lead to relations between luminosity 
and light curve shape that are opposite to the observed behavior. They conclude that the 
total mass of the explosion is a natural and simple explanation of the observations. 

If a wider range of masses could be contributing to SN la explosions, then the possibility 
arises that different stellar populations would produce different origin functions (see also 
Kenyon et al. 1993). Moreover, since stellar populations do evolve, SN la luminosities may 
depend on the particular form of the WDMF which should evolve with redshift. In principle, 
these dependencies can produce biases and selection effects in surveys for distant SNe which 
become important to evaluate and remove. These include the standard Malmquist bias 
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concerns, as well as concerns about galaxy type and position dependencies in the detected 
SN samples.^ These concerns are standard ones with any extragalactic sample, although 
they have not been extensively investigated in the SN studies to date, partly due to the 
small sample size (see Ruiz-Lapuente et al. 1995). 

At issue here is whether the progenitors of SNe la have a significant range in mass, 
that in turn produces a range in SN la luminosities, or if SNe la principally come fron 
Chandrasekhar mass white dwarfs. In this paper, we focus attention on the first possibility 
and produce a series of models which produce different white dwarf mass functions 
(WDMFs) for differing star formation histories. These models have predictive power for 
both the range in SN la luminosities as well as the mean SN la luminosity for a given mean 
stellar population age. 

However, regardless of the physics that produces SNe la, it is now well-established that 
empirical corrections to their luminosity based on the form of the light curve (e.g., Riess et 
al. 1995,1996; Hamuy et al. 1995, 1996b) produces a Hubble diagram which is linear out to 
z = 0.1 with a scatter of < 0 m .15. To first order, this argues that the intrinsic range of SN 
la luminosities is irrelevant as the multi-color light curve (MLCS) and/or luminosity-decline 
correlations empirically correct for this range. In fact, these empirical corrections may be so 
good that systematic differences between galactic stellar populations may now be revealed. 
There is already some observational evidence bearing on this (Hamuy et al. 1996) and so 
the principle task of our modelling procedure is to demonstrate how systematic differences 
in galactic stellar populations directly lead to systematic galaxy-galaxy differences in SN 
la luminosities; thus the observed galaxy correlations may have the power to discriminate 
between different progenitor models. 

In fact, we will demonstrate that our model predicts evolutionary corrections to SN la 
luminosities that, at z = 0.5, are an important percentage of the total cosmological signal 
differential between q D = 0 and q G = 0.5. In light of the concerted efforts being made in the 
detection of SNe la at redshifts > 0.3 (Perlmutter et al. 1995) and the expectation that 
fundamental cosmological parameters can be determined it seems especially important to 
understand, in as much detail as possible, the dependence of mean SN la luminosity on the 
underlying stellar population. 

The thrust of this paper evaluates the WDMF and its variation as a function of stellar 
population and evolutionary state, under the assumption that the dispersion of M^(SN 


1 For example, a radial gradient in luminosities of SNe in galaxies could easily exist because of abundance 
and age gradients, and combined with the reddening distribution, systematic luminosity differences for SNe 
in the outer regions of spiral or irregular galaxies might dominate the observed samples. 
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la) is correlated with the mass of the WD progenitor. In particular we focus on one issue: 
how does the WDMF depend on its parent stellar population (Section 2)? We then apply 
a simple parameterization of that dependence to two cases of interest: SNe la arising 
in different galaxy types (populations with different star formation histories), and the 
dependence of progenitor mass on cosmological look back time (Section 3). We evaluate 
these effects on the determination of q Q from z > 0.3 SN la detections. Our concern is in 
the scatter in the candle, and the zero point of the flux scale is irrelevant for this discussion 
(but relevant for H 0 ). 


2. A Simple Parameterization of the White Dwarf Mass Function 

We consider the luminosity distribution of SNe la to be a separable function, A, which 
can be written as 


A = G(m) N wd L(m wd ), 


( 1 ) 


where G{m ) is a source function, i.e. a restriction beyond the stellar population inputs 
on the mass range of WDs which can become SNe la and which would include various 
pathways (binary formation, etc.), N wd is the number distribution given by the WDMF, 
and L{m w d ) is the conversion from WD mass to luminosity (essentially the Ni mass core 
of the exploding WD). For this paper we will assume that G(m ) = 1 (i.e., no additional 
restrictions beyond those which we model); for the Chandrasekhar mass ignition model, 
this function would be a delta function at 1.4Mq. 

We build a simple model of the WD mass distribution as a function of population 
age based on prescriptions for the WD initial mass - final mass relation, theoretical stellar 
lifetimes, and a star formation rate (SFR) parameterization. We assume that the rate of 
formation of stars of a given mass at a given time can be characterized by a separable initial 
mass function (IMF) and SFR, as 


R(m,t) = $(m) ( A/t s ) e ^ td , (2) 

where m is mass in solar units, A is a dimensionless normalization, t s is the dimensional 
time unit, t d is the decay time of the SFR, and $(m) is the IMF of the form 


<F(m) = N 0 m a . 


(3) 



5 


The number of stars which will leave the main sequence to become WDs in a given 
mass interval, dm, and in a given time interval, dt, is 

dN evo i(m, t ) = 0(f — r(m)) R(t — r(m)) dmdt, (4) 

where r(m) is the timescale of evolution for a star of mass m, and 0 is a step function 
equal to 1 for t > r(m) and 0 otherwise, and which allows the two cases of t < r(m) and 
t > r{m) to be compactly written. To determine the number of stars which have evolved 
into WDs by a given time, t, the above equation is integrated over t to yield 

N e voi(m,t) — 0(f — r(m)) f dt$(m) ( A/t s ) e -(*- r M)/*d dm. (5) 

J r(m) 


Letting t' — t — r(m), then 


ft—r{m) 

N evo i( m ,t ) = Q(t — rim)) &(m)dm / dt\A/t s )e~ t ' td (6) 

Jo 

= 0(f — r(m)) $(m)dm A (t d /t s ) [1 — e -( t - T M)Ad]_ ^ 

To convert N evo i(m,t) to N wd (m,t) requires an initial mass - final mass relation, which 
we achieve from a quadratic parameterization of empirical relation “A” of Weidemann & 
Koester (1983): 


M wd = 0.48 - 0.016 m + 0.016 m 2 , (8) 

where m is the initial main sequence mass of a star as used above, and M wd is the mass 
of the resulting WD. This agrees with observations, which are very limited, and gives a 
1.376Mq WD for m = 8 Mq, which we assume to be the highest mass star that produces a 
WD remnant. Clearly, all results we obtain subsequently stem from this parameterization, 
and so, to the extent that it can be justified by the observations, we have a reasonably firm 
foundation. 

We also require r(m), the pre-WD evolutionary timescales as a function of mass, which 
we achieve from a re-parameterization of the equations given in Eggleton et al. (1989). We 
simplify their parameterization as we require only lifetimes for stars with masses from 1 to 
8 Mq, whereas their equations are valid for 1 to 80Mq. Additionally, we renormalize their 
stellar lifetimes so that a 1 Mq star has a main sequence lifetime of 10 10 years. Following 
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Eggleton et al. (1989) we also take the post main sequence lifetime to be 15% of the main 
sequence lifetime. The resulting parameterization is then 

r(m) = t Q m~ 2 ' 8 , (9) 

with t Q = 1.15 x 10 10 years, for m — 1 - 8Mq. For an 8Mq star, the timescale of evolution 
is 34 Myrs. 

To use the above equations we set A = 1 (arbitrary normalization) and t s — 1 Gyr (i.e., 
all time units in Gyrs). We then choose various SFR models with td = 1,3,5,10, and 100 
Gyrs to simulate the range from single age ellipticals to constant star formation spirals. We 
explore the range a = 0 to — 3 (a = —2.35 is the Salpeter value) for the slope of the IMF 
and then calculate N evo i over the range of time until t = 12 Gyrs. Finally, we transform 
N evo i to N w d via the initial mass - final mass relation. 

Figure 1 shows the resultant WDMF for different values of the mass function slope 
(i a = —3, —2, —1,0) in each panel, for two different SFR decay times (1 Gyr, essentially a 
burst; and 100 Gyrs, almost constant SF), and for 6 different ages (0.5,1, 2,4, 8, and 12 
Gyrs) since the onset of SF. 

One immediate test of these models is a comparison with the solar neighborhood WD 
mass function. Figure 2 shows the observed mass function of Bergeron et al. (1992). As 
they note, this magnitude-limited survey selects against the fainter, low radius (high mass) 
WDs. Additional selection may also be caused by the quicker cooling of higher mass WDs, 
and possible scale height inflation that would preferentially select against all stars of higher 
mass than the current turn-off mass of the disk population. Plotted on this distribution is 
our 10 Gyr, steady SFR, a = —2.35, model with an arbitrary normalization. The agreement 
is satisfactory after noting that Bergeron et al. (1992) interpret the lowest mass WDs (first 
several bins) as likely results of binary evolution. On this basis, we believe that our models 
produce WDMFs that are astrophysically reasonable. 


3. Luminosity Functions: Predictions, Samples, and Biases 

Woosley & Weaver (1994) explored the details of 0.6 - 0.9Mq WDs accreting from 
a companion post main sequence star. They found a number of scenarios where 0.1 to 
0.2 Mq of material (He) could be accreted before a thermal runaway in the surface layers 
occurred. Since these thermal runaways propagate more rapidly around the surface of the 
WD than the resulting shock wave propagates into the interior of the WD, the shock wave 
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is focused in the deep interior, often resulting in a detonation. We use their models as 
the basis of our parameterization of the amount of light given off by the supernova (based 
on 56 Ni production) as a function of mass of the accreting white dwarf. The Woosley & 
Weaver models are meant to explore a range of accretion rates and metallicities, and we 
parameterize their results as model A, which has a mass accretion rate of 2.5 x 10~ 8 Mq 
yr~ x , and model B, which has a mass accretion rate of 3.5 x 10 _8 Mq yr _1 . Model A creates 
SN la type explosions for a pre-accretion mass as low as 0.6Mq, whereas Model B creates 
SNe la for masses as low as 0.7 Mq. The upper mass limit of their pre-accretion WDs is 
0.9 Mq, but we will assume that this relation can be extrapolated up to 1.1 Mq, which is 
a likely upper limit to C-0 WDs (Iben & Webbink 1989; but see Kippenhahn & Weigert 
1990). Our extrapolation and the unknown upper limit is overly simplistic, but is sufficient 
for our purposes. Our parameterizations of these two models are then 


L a (m wd ) ~ m™{m wd ) = -1.2 + 2.4 m wd for m wd > 0.6 and (10) 

Lb{m wd ) ~ m™(m wd ) = -1.3 + 2.3 m wd for rn wd > 0.7, (11) 

where 56 Ni masses in excess of 1.376Mq are set to 1.376Mq. This adjustment only affects 
WDs in the incremental mass range 1.07 - 1.10Mq, and only for model A. 

The resultant luminosity functions, A, from the product of N(m ) and L(m wd ), are 
shown in Figure 3 for several combinations of age, a, and SFR parameterizations. The A 
functions are very flat, as expected from the nature of the almost power law mass functions 
and the simple linear relation between the Ni mass and the WD progenitor mass. These 
luminosity functions are strongly non-gaussian, which is likely to be the result in general if 
the wide mass range assumption we have made (essentially the Woosley & Weaver models) 
are not given any features by the source function, G(m). 

Ideally we would now like to rigorously compare Figure 3 with the observed SN la 
luminosity function. However, it is our contention that the observed LF is poorly known. 
Essentially all extant surveys have to be corrected for completeness. These incompleteness 
corrections depend on the assumed intrinsic form for APj. As such these corrections usually 
have the flavor of self-fulfilling prophecies: the derived “a” will depend on the assumed 
dispersion. Since the number of well-measured SNe la occurring in host galaxies with well 


2 So A obs = 5(A), where 5 is a selection function which depends on the magnitude limit and other 
properties of the survey. 



measured distances is quite low (e.g., the 9 objects in Phillips 1993), neither the intrinsic 
LF nor a reliable estimate of the mean M^(SN la) can be made from extant data. 

As a result of data paucity, the construction of the proper SN la LF is currently 
an ambiguous and contentious issue which remains unresolved. Figure 4a shows the B 
luminosity function for 29 SNe la from the Calan/Tololo survey (Hamuy et al. 1995, 1996a), 
plus the 9 objects from Phillips (1993). One of our referees argued that the 9 SNe from 
Phillips (1993) “over represents” peculiar SNe la. We feel this reasoning was circular since 
the criteria for inclusion in the Phillips sample is only that a good distance to the galaxy 
has been derived (from surface brightness fluctuations or Tully-Fisher measurements). 
How could selection based on the existence of an independent distance estimate cause an 
over representation of anomalous SNe la? In fact, the Phillips criteria is exactly what 
should be used in the construction of a representative LF as long as no identifiable bias 
exists in the distance determinations to these 9 galaxies. We also choose the Calan/Tololo 
sample because it is the largest collection of SNe la with homogeneous (although still not 
quantified) selection criteria. Figure 4a presents the observed LF, uncorrected for any 
probable selection effects. The Phillips (1993) sample of nearby SNe la in galaxies has 
unknown selection effects, while the Calan/Tololo sample of southern SNe has some galaxy 
type dependencies with distance that are still being explored. For this sample, absolute 
magnitudes of SNe la are assigned using redshift as the distance indicator. The most 
significant aspect of Figure 4a is not its shape or mean value but rather the total luminosity 
range that is exhibited. The very faint object evident in this figure is SN1991bg, a very red 
SN that has been universally tagged as being an anomalous SN. 

Figure 4b shows the SN la LF for all the SNe from Vaughan et al. (1995) (hereinafter 
VBMP) with data obtained after 1970 (30 SNe). If we exclude from the first sample 
SN1991bg, the two distributions are essentially identical. For the 30 objects in Figure 
4b, the mean B-magnitude is —18.50 ± 0.49 while the mean B-magnitude for the 37 
objects in Figure 4a is —18.50 ± 0.4. These dispersions are relatively large, and obviously 
uncorrected SNe la would not appear to be a premier distance indicator. VBMP claim to 
be able to lower this dispersion by identifying and removing SNe with deviant red or blue 
color. Since the intrinsic spectral energy distribution of SNe la is not yet well known from 
theory, a color-based rejection criteria is at best risky. If we examine the 20 objects in 
the VBMP sample that were discovered after 1980 and reject SN1991bg and SN1986G (as 
obvious deviants), the mean magnitude is —18.40 ± 0.46 (a change of —0. m l in the mean 
is significant in the cosmological context). VBMP reject 3 more objects from this sample, 
including the very well studied object SN1989B. VBMP specify the observed color (B—V = 
0.30) of SN1989B as being anomalous but Wells et al. (1994) attribute its color to a large 
reddening, specifically E(B—V) = 0.37. Removing this single object from the 18 most recent 
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SNe in the VBMP sample lowers the dispersion from 0.46 to 0.34 mag! Yet if reddening 
is the reason for the anomalous color, then obviously the absolute magnitude of SN1989B 
is substantially brighter than the value listed in VBMP. After trimming of the anomalous 
objects in Figure 4b, VBMP find a distribution with a mean magnitude of —18.54 ± 0.35. 
This mean is very similar to the values we derive for Figure 4a. 

We contend that the SN la LF is simply not yet well-determined due to limited 
sample sizes and survey volumes, and the difficulty of determining direct and independent 
distances to many of the host galaxies. While we may have a reasonable estimate for 
the maximum brightness of SNe la, we do not know the entire LF. Furthermore, the SN 
la LF of VBMP is not representative of the whole distribution of SN la luminosities but 
rather of VBMP’s selected sub-sample in which they have chosen to ignore or reject a 
fair percentage of the faintest SNe. Are these rejected objects not, therefore, SNe la? 
Without an adequate explanation of the mechanism that causes the rejected objects to be 
anomalously underluminous, it seems premature to a priori exclude them when specifying 
the intrinsic range of SN la luminosities and then claim that the sample of distant SNe la 
is identical to the selected sub-sample. 

For example, we know that the Calan/Tololo sample has obvious selection effects; SNe 
la fainter than —18.0 will not be found in at least half the surveyed volume, since they 
fall below the apparent magnitude cutoff of the survey. Some of these selection effects are 
discussed in Hamuy et al. (1994). While we are not prepared here to analyze completeness 
of extant SNe samples, we schematically illustrate our concerns in Figure 5, which shows 
the redshift distribution of the Calan/Tololo SN la sample. This distribution is extremely 
flat, with a median recessional velocity of ~ 14,000 km s' 1 , but extending out past 30,000 
km s' 1 . We have included lines in this figure to demonstrate the expected increase in the 
sample due to volume effects. The dashed line is a normalization assuming the survey is 
complete out to 4,000 km s' 1 (which may be representative of the lower luminosity SNe 
la), while the dotted line assumes it is complete out to 14,000 km s' 1 (representative of the 
brighter SNe la). In either case, we conclude that the sample is severely incomplete through 
much of its volume; in the first case it is 98% incomplete at the median redshift of 14,000 
km s' 1 . If the deeper completeness normalization is assumed, the excess of low redshift 
SNe are those of fainter absolute magnitude. These comments about incompleteness in the 
samples are a reflection of our concerns about the completeness and accuracy of the extant 
SN samples which have been used to construct the SN la LF. 

Our simple model of the range of SN la luminosities attempts to explore the systematic 
connection between this range and galaxy type. Indeed, it may be very difficult to judge if 
the very distant SNe la have a broad or narrow distribution of absolute magnitudes, due 



to selection effects and cosmological corrections. We are exploring a scenario that makes 
certain predictions that can be tested on both local and distant samples. 


3.1. Galaxy Population Dependencies 

For the purposes of this discussion, we consider that our simple models may be 
assigned to galaxy types based on stellar population type. Thus we will assume that we 
can characterize elliptical galaxies as single burst models, with the majority of stars formed 
12 Gyrs in the past, and with a 1 Gyr exponential decay time. We will also assume a IMF 
slope of a = —2 for this single burst model. We will assume an actively star forming galaxy 
(SFG) can be characterized by star formation starting approximately 8 Gyrs ago (the 
approximate age of the Galactic disk) with effectively continuous star formation (G = 100 
Gyrs) and with several IMF slopes. Clearly, these assumptions can be challenged, but they 
correctly predict the average UBV color differences between spiral and elliptical galaxies 
(see Larson & Tinsley 1978; Bothun 1982). For the cosmological parameters we will assume 
H a = 50 km s -1 Mpc -1 (after all this is a paper on SNe) and q Q = 0.5. This universe is less 
than 14 Gyrs old. 

Although the power law nature of the LFs shown in Figure 3 precludes the calculation 
of a reliable mean SN la luminosity as a function of galaxy type, we can use these means 
to make rough estimates. Figure 6 shows the mean SN la luminosities that are obtained by 
integrating over the Ni mass distribution as normalized by total number of SN la events. 
This figure demonstrates that our model LFs do not change shape after ~ 0.5 Gyrs, which 
is the evolutionary timescale of the lowest mass progenitors (3Mq) which explode as SNe 
la. We caution, however, that our model does not include the unknown, but possibly large, 
time delays inherent in the binary mechanism before mass transfer begins. Thus the large 
change in the SN la LF shape evident at early times should take place over a greater time 
period, making the SN la LF more sensitive to stellar population age than this figure 
implies. Figure 6 demonstrates that, in the case in which ellipticals and spirals have the 
same IMF slope a, the differences in mean SN la luminosity are small (< 0.04 mag). This 
small difference is not surprising as star formation histories with similar a will produce 
very similar WDMFs once the mean age is greater than the evolutionary timescale of the 
WD progenitors. The difference in mean SN la magnitude between an a = —2 elliptical 
and an a — 0 spiral is 0 m .22. (An even more extreme difference occurs with large lookback 
times.) The data from the SN samples support these dependencies on the underlying stellar 
population. Hamuy et al. (1996a) find that the brightest SNe occur in late type galaxies 
(see their Figure 3) and even more strikingly, that a strong correlation exists between the 
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decline rate of the SN la light curves and the host galaxy morphology (see their Figure 
4). Branch et al. (1996) find that SNe la that occur in “red” galaxies are 0.3 magnitudes 
less luminous than those that occur in “blue” galaxies. While our model calculations are 
meant to be illustrative only, they do show that differences in mean stellar population age 
and/or slope of the IMF can produce significant differences in mean Ms(SN la) that are 
approximately the same size as the effects seen in existing data samples. 

However, variations in mean M#(SN la) between galaxy types are not the relevant 
quantity with respect to distance measurements. Rather, the total range of M^(SN la) is 
important, especially when considering the effects of Malmquist bias. In Figure 7 we plot 
the initial mass - final mass relation (equation 8) together with nickel mass production. The 
models indicate a range of 4 - 5.5 in Ni mass, which indicates a range in SN la luminosities 
of up to l m .8. We terminate our masses at a 1.4Mq WD, but if WD mergers at all masses 
are a possible channel, then the functions should be continued up to the possible sum of 
2.8 Mq, giving a total possible range of 2.6 mags for SN la luminosities. 

Our predicted range is similar to the observed range shown by Phillips (1993). However, 
this comparison is only indirect. Our predicted results are for the SN la luminosity range in 
a single galaxy for a specific WDMF, whereas in comparing to observations, we are sampling 
over a range of galaxy types. Still, the rough agreement between the results based on our 
model parameterization and the available observations has an alarming implication: an 
order of magnitude range in M#(SN la) immediately suggests that Malmquist corrections 
are large for any distant extragalactic sample of SNe. We believe the current observations 
are effectively sampling this range. That the method of Riess et al. (1996) can lead to 
Hubble diagrams with such low dispersions indicates that the light-curve corrections to SN 
la luminosities are very effective at compressing this intrinsic luminosity range. If these 
corrections continue to work well in larger samples, then it becomes clear that the intrinsic 
luminosity range of SNe la is essentially irrelevant with respect to determining distances. 
All that is required is a secure calibration of these light-curve corrected luminosities. 

For the simple case in which the Ni mass is proportional to the WD mass (e.g., 
equations 10 and 11), our models predict a spread in SN la luminosities of at least l m .5. 
However, invoking variations in WD progenitor mass as the sole cause of the SN la 
luminosity spread may not be necessary. Various explosion scenarios can easily give 50% 
variation in the energy release (Khokhlov et al. 1993; Hoflich et al. 1996). Whether these 
models would have a systematic dependency on the WDMF or on the evolutionary state 
or metal abundance undoubtedly depends in detail on the nature of the explosion. We 
have also ignored any effects of changes in the WDMF on the binary frequency or whatever 
progenitors are the SN la source G(m ) (cf. Kenyon et al. 1993). While the wide binary 
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source function may be independent to zeroth order of the details of the individual stellar 
mass function, the SN la source function probably evolved in a more complicated manner 
than a simple dependence on mean WD mass. Hence, several physical effects can cause the 
SN la LF to depart significantly from a delta function. 


3.2. Cosmology and Evolutionary Corrections 

The redshift - magnitude relation in standard form yields the equation 

m = M + 25 — 5 log H 0 + 5 log cz + 1.086 (1 — q Q ) z + .... (12) 

At z= 0.5, for H a = 50 km s^ 1 Mpc -1 , the difference between q Q of 0.0 and 0.5 (empty 
versus critical mass models) is 0. m 27 (assuming a zero cosmological constant). With 
photometric accuracies of ~ 0. m l per SN event, statistics of a sample of 10 well-measured 
objects would permit an ~ 8cr discrimination between empty and critical models. Clearly, 
however, systematic errors or biases at the 10% level become very significant and lead to an 
effective q a measurement .[] 

For these cosmological parameters the look-back time at z = 0.5 is 3.75 Gyrs. 
Inspecting our models we see little change in our E galaxy sources. However, the actively 
SFG shows significant evolutionary effects in the sense that the WDMF is populated toward 
the more massive objects in the past, and thus the A function produces brighter SNe. The 
WDMF also depends sensitively on the assumed age parameter for the SFG; if the galaxies 
are assumed to initiate star formation 8 Gyrs ago, at z = 0.5 the mean luminosity can be 
as large as 0 m .31 brighter for the a = 0 case.[] 

The nature of the general galaxy population at z — 0.5 is still poorly determined. At 
a minimum, however, we expect the fraction of young or starbursting galaxies in some 
random field to be significantly higher than is observed at the present epoch. The SN la 
rate at some epoch, z, is a function of the total number of young stars that exist at that 
epoch since massive white dwarfs are produced by short-lived stars. Hence, a single, massive 


3 We would measure q e J t = q a — (dL/dt)/L/H a , in the case of luminosity evolution, for example. 

4 We ignore here the K-correction issue, which can be complicated at the few percent level for complex 
spectral types such as SNe. For example, different K-corrections are probably necessary at different epochs 
for a given supernova event, because of the significant evolution of the spectral energy distribution. Clearly, 
shifting the observed bandpass with redshift is an important consideration (cf. Hamuy et al. 1993; Kim et 
al. 1996). 
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starbursting galaxy could completely dominate the rate. Additionally, early star formation 
may be characterized by a far different IMF slope than we observe today. This would be 
particularly troublesome as the SN la LF changes dramatically with IMF slope (see Figure 
6). For our purposes, we pose one specific question: how much of a star formation burst 
is required to create a significant number of SNe la from the burst population relative 
to a single-age 5-8 Gyr (elliptical) galaxy? SNe that occur in ellipticals (or the old 
population in a spiral which we assume is negligible at these redshifts) can be thought 
of as being the background SN population against which SNe occurring in star bursting 
galaxies are detected. The relative contribution of the SF and background populations can 
be parameterized as 

No. burst SNe _ n burst exp{-t/t 0 ) _ n burst exp(-l or - 0.5) 

No. background SNe n back exp{—t/t 0 ) n back exp( —8 or — 5) 

For population age, t = 5 Gyrs, 

(SN burst /SN back ) {o burS fjn back ) (55 90) (14) 

and for population age, t = 8 Gyrs, 

( SN burst /SN back ) = (n burst /n back ) (1100 - 1800). (15) 

Table 1 provides the percentage of SNe resulting from a burst as a function of the 
percentage of mass in the burst population relative to underlying stellar population. Again, 
the values in this table come from considering one single age elliptical with one starbursting 
spiral. Column 1 lists the percentage of total SNe that come from the starbursting spiral, 
while columns 2 and 3 list the starbursting mass fraction for the t = 5 and 8 Gyrs cases. 

For the case of low burst strength (e.g., < 2%), we find the expected result that since 
the field contains two galaxies, 50% of the SNe come from one of the two galaxies. However, 
for a star formation burst of 10 - 20%, 90% of the SNe will come from that one starburst 
galaxy. The situation is even more extreme if we consider a true starburst galaxy (burst 
strength > 100%) in which case 99% of the SNe come from that one galaxy. These results 
indicate that if a field at z = 0.5 contained 90% ellipticals and 10% starburst spirals with 
burst amplitudes of 10 - 20%, then 50% of the total SNe generated by these galaxies would 
come from the minority population. If, however, these galaxies are preferentially dusty, 
then extinction effects may reduce the detection of SNe la from these hosts. Thus, samples 
of distant SNe might be dominated by host galaxies which have young mean ages. 



-14- 

4. Conclusions and Caveats 

Overall, our simple model of the dependence of the SN la luminosity on the underlying 
WDMF allows us to make the following predictions: 

1. In the mean, the SNe la occurring in spiral galaxies should be more luminous than 
those occurring in elliptical galaxies; bright SNe la in E galaxies should be very rare. This 
effect can be seen in the data compilations of Phillips (1993) and in the Calan/Tololo survey 
(Hamuy et al. 1995, 1996a). 

2. A correlation should exist between SN la luminosity and the color of the host galaxy 
population, with brighter SNe present in bluer galaxies. 

3. The SNe in the disks of spiral galaxies should be more luminous in the mean than 
those in the bulges. This effect may be hard to observe because of reddening effects. A 
reddening independent light curve parameter (such as Amis) should correlate with position 
in a spiral galaxy, with the broader light curves (smaller Am 15 values) preferentially in the 
disks or outer regions of the spirals. 

4. More distant SNe should show slower light curve decay (smaller Amis values) than 
the nearby sample because these SNe are preferentially more luminous. This prediction is a 
consequence of both starbursting galaxies dominating distant samples and the Malmquist 
bias that directly results from the large range in intrinsic SN la luminosities. 

5. Distant SNe are expected to come predominantly from bright, blue, spiral or 
irregular galaxy hosts, most of which are in an elevated state of star formation. The mean 
age of these hosts will be younger than the mean age of most z = 0 calibrating galaxies, 
making it important that starburst galaxies like NGC 5253 are included in the local 
calibrating sample. We have already shown that M B (max) is sensitive to the mean age of 
the stellar population. Thus correcting for this mean age effect requires detailed knowledge 
of the nature of the stellar populations in distant galaxies. The predicted difference in 
Ms (max) obtained under modest assumptions about the star formation history of galaxies 
is an appreciable fraction of the cosmological signal that distinguishes q 0 — 0 from q Q = 0.5. 

6. The form of the LF for SNe la should be approximately a power law (see Figure 3), 
if we assume the binary formation function introduces no strong features. Larger samples 
of low redshift SNe will be needed to determine this function. 

7. In general, the form of the WDMF predicts a range of 1.5-2.5 m in SN la luminosities. 
We have argued that current SN la samples have effectively sampled this range in luminosity 
and that their usefulness as a distance indicator depends critically on the universality of 
the light curve correction algorithms (e.g., MLCS) in compressing this luminosity range. 
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If many of these predictions are borne out, we would contend that such observational 
evidence favors the sub-Chandrasekhar mass hypothesis as the main SN la progenitor. In 
fairness, our results and modelling procedure and its application to the SN la distance scale 
are subject to a number of caveats and we close this paper by discussing them. 

In converting WD masses into SN la luminosities we use the recent calculations for 
sub-Chandrasekhar explosions by Woosley & Weaver (1994). These models accurately 
reproduce the observed correlation between decline rates of the light curves and luminosity, 
and are able to produce more 44 Ti and 48 Cr than other types of models (see discussion in 
Livne & Arnett 1995), which is important in matching solar abundances. The ID treatment 
of Woosley & Weaver (1994) yields similar results to the 2D treatment of Livne & Arnett 
(1995). We choose the Woosley & Weaver models because they have clear predictive power, 
not because we consider these models to be definitive. While models of Chandrasekhar mass 
explosions can also yield a range of luminosity (e.g., Hoflich et al. 1995; 1996), based on the 
nature and degree of turbulence in the explosion, we contend that, because the dispersion in 
SN la luminosities is not small and seems to be correlated with galaxy morphology, effects 
in addition to explosion physics most likely produce the observed LF. It is not our intention 
to delve into SN explosion physics or discuss which SN models in the literature are more 
nearly correct. Instead we have argued that a major part of the observed luminosity range 
for SNe la can result from a dependence of mean SN la luminosity upon the mean stellar 
population of the host galaxy. 

In this case, the mixture of host galaxy types in any SN la sample determines the 
LF for that sample. Thus, it is not surprising that there is disagreement over the form of 
a typical host galaxy in the SN la sample at z = 0. Moreover, our models clearly show 
the importance of starbursting galaxies in distant samples. The higher SN la rate in these 
galaxies allows the minority population to dominate the observed frequency. Since these 
galaxies have younger mean ages and hence more extended WDMFs, the range of SN la 
luminosities is larger than that in a z — 0 spiral or elliptical. Indeed, the distribution of 
SN la host galaxies in the nearby Universe shows some curious properties which makes it 
hard to determine if the typical host is a spiral or an elliptical. For instance, the modestly 
star forming galaxy M 100 has had four detected SNe since 1901 (one of which is a type 
II), whereas the megastar elliptical M 87 has had zero. NGC 5253, a low mass but actively 
star forming galaxy, has had two detected SNe la in the last 100 years. By comparison, the 
Coma cluster, home to ~ 10,00 gas poor L* galaxies (e.g., 10 4 NGC 5253 masses), has not 
had a single SN la event detected for the last 22 years. 

In contrast to this anecdotal evidence, which suggests that spiral hosts dominate over 
elliptical hosts, the 30 or so SNe la that have been detected in the Calan/Tololo survey 
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show nearly equal numbers of elliptical and spiral hosts beyond z = 0.033, demonstrating an 
anff-Malmquist bias. The dominance of nearby spiral hosts [at redshifts z < 0.033 (Hamuy 
et al. 1996a)], may be a result of the avoidance of nearby clusters in the search fields. 
However, in the distant half of the sample no a priori selection against clusters existed, 
and hence, proportionately more ellipticals should be in that sample, causing some of the 
variation. Thus variation in the S/E host ratio could reflect these selection criteria, as well 
as the low space density of relatively unreddened starburst spirals in the local universe 
(z < 0.1). ft is unlikely, however, that similar circumstances would continue to hold at 
larger redshifts. 

Finally, we comment on the use of light curve corrections to SN la luminosities in the 
context of our model. Astrophysical measurements, based on either light curve parameters 
(Phillips 1993; Riess et al. 1995,1996; Hamuy et al. 1996a) or spectroscopic analysis (Nugent 
et al. 1995), appear to correlate well with peak SN luminosity. Hamuy et al. (1996a) show 
that high quality light curves exhibit a characteristic shape and luminosity - decay time 
relation (Phillips 1993; hereafter the Phillips relation) that produce significantly improved 
peak magnitudes and much more accurate relative distances than the use of a single 
absolute magnitude calibration (see also Riess et al. 1995,1996). In the model explored 
here, the Phillips relation represents a stellar mass sequence. Although we do not attempt 
to derive the relation between light curve parameters and mass explicitly, the luminosity - 
mass relation is itself linear. 

Intrinsic dispersion around the Phillips relation or the MLCS relation of Riess et al. 
(1996) would be caused by additional parameters (e.g., metallicity), which may or may not 
correlate with the host galaxy stellar population and the WDMF. It is the dispersion around 
these relations that become directly relevant to correcting distant samples for Malmquist 
bias. Moreover, corrections to SN la peak luminosities using z — 0 light curves may not 
be strictly applicable to distant SNe in star-bursting galaxies as such galaxies will be rare 
(or perhaps non-existent) in the nearby calibrating sample. At the very least, our models 
show that knowledge of the WDMF as inferred from the nature of the stellar population 
of the host galaxy is critical in order to determine potential systematic differences in 
light-curve shape between the distant host galaxy and the calibrating sample. Minimizing 
these differences may well validate the approaches of Hamuy et al. (1995, 1996b) and Riess 
et al. (1995,1996) that have produced linear Hubble relations out to z — 0.1 with a scatter 
of approximately 0 m .13 - 0 m .17 (Hamuy et al. 1996b, equations 7-9) which raise the 
expectation that q Q can be determined from such data. 


We acknowledge helpful discussions with Karl Fisher, Philip Pinto, and Michael 
Richmond. We also acknowledge Mark Phillips, Mario Hamuy, and Nick Suntzeff for 



inspiring us to investigate the possible connection between SN la luminosities and the 
underlying stellar population. Finally, we wish to dedicate this paper to Marc Aaronson, 
who would have wanted a thorough investigation of the reliability of SNe la as standard 
candles. 
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Fig. 1.— WDMFs for different IMF slopes (panels display exponential SF duration and IMF 
slope) and for 6 different ages (0.5,1,2,4, 8, and 12 Gyrs, bottom to top, respectively), a) 
WDMF of burst models with 1 Gyr exponentially decaying star formation; b) WDMF of 
active star forming galaxies with 100 Gyr exponential star formation. 

Fig. 2.— Observed Galactic WDMF (from Bergeron et al. 1992) plotted versus 10 Gyr, 
a = —2.35 model for an active star forming galaxy (100 Gyr exponential star formation). 

Fig. 3.— The differential SN la luminosity functions of the cumulative (time integrated) 
stellar population as derived in the text, for the burst models (a, b) and active SFGs (c, d) 
for 2 different mass function slopes (a = —2,0) and 5 different ages (1,2,4, 8, and 12 Gyrs, 
top to bottom, respectively). 

Fig. 4.— a) Histogram of absolute B magnitudes for SNe la in the sample of Phillips (1993) 
and the Calan/Tololo survey (Hamuy et al. 1996a). The magnitudes of the 9 events from the 
Phillips paper have been adjusted to a zero point consistent with the Calan/Tololo sample 
calibration, b) Histogram of absolute B magnitudes for the SNe la after 1970 from the 
sample from Vaughan et al. (1995). 

Fig. 5.— The velocity histogram of the Calan/Tololo SN la sample. Shown for comparison 
are lines showing the volume increase, normalized to the observed SN number at cz = 4,000 
km s' 1 (dashed) or cz = 14,000 km s^ 1 (dotted). Under either assumption, the sample is 
seriously incomplete (see text). 

Fig. 6.— The mean SN la luminosity as a function of age, SF history, and IMF slope for 
our model A. The luminosity is in terms of equivalent Ni mass. The clustered lines at a 
particular IMF slope are for models with different SF histories, with the solid line being the 
t d — 1 Gyr (burst) model, the dashed line being the t d = 10 Gyr model, and the dotted line 
being the t d = 100 Gyr (steady SF) model. 

Fig. 7.— The relevant mass ranges for creating model SNe la. The line marked “WD’ : 
shows the initial mass - final mass relation, whereas the other two lines display the Ni mass 
produced in the SN la explosions for models A and B, as a function of the main sequence 
mass of the progenitor. 
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Table 1: Burst statistics 


N from burst 

ft burst/^back t 5 

ftfeurst /ft-bacfc ftt t 8 

50% 

1- 2% 

< 0.1% 

90 

10- 16 

< 1 

99 

110-180 

6 - 9 






























